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SUMMARY 

A statistical, theory "based on the presence of small crystalli- 
zation nuclei suspended in water is developed to explain experimental 
results showing that on the average small droplets can he supercooled 
to lower temperatures than large ones . Small nuclei of crystalliza- 
tion are assumed responsible for causing supercooled water to freeze 
spontaneously. 

The average behavior of supercooled droplets is reproduced on the 
basis of probability theory with an assumed distribution of crystalli- 
zation nuclei with respect to the temperatures at which the nuclei 
cause freezing. The most probable distribution curves of spontaneous 
freezing temperatures for water droplets of various sizes within the 
size range found in clouds are obtained. 


INTBODUCTION 

The prolonged operation of unprotected aircraft in supercooled 
clouds has been found to be a hazard. For the icing hazard, to be over- 
come, the physical conditions under which supercooled clouds exist must 
be understood. The physics of supercooling" and freezing is incompletely 
understood, as shown by conflicting opinions expressed in the litera- 
ture. The following factors have been described as causes of freezing 
in supercooled water samples: (1) introduction of ice crystal; (2) mech- 

anical shock, such as vibration, scratching of solid surfaces together 
in water sample, and, sound waves; (3) length of time water has existed 
in supercooled state; (4) impurities; and (5) droplet size. 

Investigators seem to agree that an ice crystal introduced into a 
sample of supercooled water will oause the water to freeze. Eecent 
experiments by Dorsey (reference 1) indicate that mechanical shock, 
however, is of little importance in freezing of supercooled water; like- 
wise, time apparently is not a faotor. From the extensive investigation 
of reference 1, it is concluded that smal .1 solid particles suspended in 
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the water or irregularities on the container surface cause the super- 
cooled water to freeze at a given temperature. These hypothetical par- 
ticles are called motes and the temperature at which a sample freezes 
without any apparent - external stimuli is called the spontaneous freezing 
temperature . The spontaneous freezing temperature was found In most 
cases to he a characteristic of each sample, Which led to the mote 
hypothesis. 

The observations described in reference 1 have been substantiated 
somewhat in reference 2 in which a statistical investigation of water- 
droplet size as a factor In determining the spontaneous freezing tem- 
peratures of the droplets is presented. By observation of a large number 
of droplets, the following results were obtained: (l) The frequency- 

distribution curve with respect - to spontaneous freezing temperature of a 
given droplet size was found to have a marked peak or mode; and (2) the 
average spontaneous freezing temperature of eaoh size group decreased 
with decreasing droplet size. 

In reference 3, the conclusions are made from observation of a 
number of droplets that the_ spontaneous freezing temperature of droplets 
from 400 microns down to approximately 50 microns rapidly decreased with 
size and that droplets larger than 400 microns had a practically constant- 
freezing temperature . The results of reference 3 are qualitatively 
verified in reference 2 but only to the extent that the average freezi ng 
temperature for a given size was found to be a function of the size. 

Similar results of a qualitative nature have been obtained with 
substances other than water. In reference 4, the possible degree of 
supercooling In sulphur droplets was found to increase with decreasing 
droplet diameter. In addition, it was found more difficult to super- 
cool sulphur droplets on a metal surfaoe, such as brass or alumi num, 
than on a glass surface. The effect of droplet size on spontaneous 
freezing temperature was ascribed to the fact - that smaller droplets are 
less susceptible to external Influences than larger droplets. The sur- 
face was thought to have an influence on the crystallizing forces within 
the supercooled droplet. Substances such as gold, platinum, rhodium, 
mercury, and other metals have been supercooled to a greater extent In 
the form of small particles than in bulk (reference 5). The effect of 
droplet size on spontaneous freezing temperature is ascribed in refer- 
ence 5 to the presence of small particles of impurity In the liquid. 

When the liquid is broken up into small droplets the foreign particles 
are Isolated in a few drops allowing the remainder to supercool to the 
temperature at which molecules of the liquid combine to form crystals 
without the aid of motes. 
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The hypothesis of reference 5 combined with the mote concept of 
reference 1 indicates a -way of explaining the statistical results pre- 
sented in reference 2. A statistical theory that accounts for the two 
conclusions of reference 2 on the basis of motes or any other agent 
causing nuclei of crystal growth to form was developed at the HACA Lewis 
laboratory and is presented herein. 


THEORY 

An assumption is made that a large number of motes are present in 
liquid water. Also, each mote is assumed to be associated with a defi- 
nite spontaneous freezing temperature . The freezing temperature of a 
water sample is governed by the mote in the sample that is associated 
with the highest freezing temperature . 

In the experiment reported in reference 2, water was condensed 
from the atmosphere in the form of frost on a polished metal surface 
by lowering its temperature to approximately -35° C. The frost was sub- 
sequently melted and droplets formed on the surfaoe . Repeated applica- 
tion of this procedure made a large number of droplets of various sizes 
available for determining spontaneous freezing temperatures. Because 
the water condensed from the atmosphere was a very Bmall fraction of 
the total amount of water available in the form of vapor, the process 
was practically the same as drawing small samples of volume AT from 
a large reservoir of liquid water of volume T. Because motes of all 
kinds are assumed to be present everywhere, it is reasonable to con- 
sider them to be present in the same relative proportions in the hypo- 
thetical volume of liquid water T. 

The problem to be solved is the determination of the most probable 
frequency -distribution curve of spontaneous freezing temperatures sub- 
ject to the following conditions: 

(1) There are Nq, H^, Rg* ^3> • • • • * • mo'tes of kinds 

0,1,2, 3, ...T, respectively, present in the volume T, 

where T is the absolute value of the temperature in °C at which 
the mote of kind T causes water to freeze. 

(2) The smal l volumes that are drawn from V are of volume AT, 
and the total number of volumes AT obtainable from T is T/AT = p. 

(3) The number of draws in each set of draws from T is so small 
that p is practically unchanged while a set is being drawn. 



4 


NACA TN 2234 


First the motes of kind T are considered. There are Nrp motes 

available -for distribution among the y small volumes AV, into -which Y 
may be divided. The total number of ways of distributing Kj objects 

among y cells is Each way of distributing objects may be 

called a complexion. Next, regarding the group of all possible com- 

N +1 

plexions, there are y cells in each complexion and therefore y T 
cells in all the complexions. The probability of drawing a small volume 


in one draw with no motes of kind T in it is ■— ■ v > where Q is the 

Nrp+1 

Nm+1 ^ 

number of cells that among the y™ cells in all the possible complex- 
ions contain no motes. The probability of drawing a small volume contain- 

Q 

Ing one or more motes is 1 - " mJTT]"* 


A simple example should serve to make the foregoing logic dearer. 
Two motes are supposed to exist in a volume of water Y . From Y a 

volume AY = — V is to be drawn. The probabilities that in a single 

draw one or more motes will be in the half volume drawn or that no motes 
will be contained in it are to be determined. For solution of the prob- 
lem, the volume Y may be schematically represented by two squares of 
equal area side by side . Thus, if the moteB are numbered 1 and 2 
for identification, the complexions may be pictured as follows: 



[2~TT 


m 


1*2 


The number of possible complexions is = 2^ = 4 and the total number 

of volumes or cells in all the complexions is y =*2 =8. The num- 
ber of cells containing no motes is two and the number containing one or 
more motes is six. The probability of drawing a cell with no motes in it 


is therefore 


Q 


, n T +1 



and with one or more motes in it. 


1 


Q 

Nm+1 
h X 


3 

I * 


In the preceding example a particular case has been evaluated. 
General expressions for the probabilities must now be obtained in terms 
of ,y and Kp. As in the previous example, oells are lined up side by 
side numbered from 1 to y. In a given complexion there are Kji i, 

H T,2» • • • K T,y motes in the cells 1, 2, . . . y, respectively, where 
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the second subscript is the cell number and Hip is the number of motes 
in the i - *** 1 cell; 


1 

2 

3 

4 

5 

i 

N T.l 

W T,2 

%,3 

%,4 

%,5 


%,-i 



The total number of motes causing freezing at a temperature of -T° C is 


% = N T,1 + N T,2 + • » • + %,p. 

The total number of complexions having Nj ^ equal to 0 is (p.-l) , 

bedause Hp objects £Lre restricted to distribution among p.-l cells 


N. 


by the condition ^=0 . Likewise, there are (p.— 1) ^ complexions 
subject to each of the conditions 2 =0, . . . Hp ^=0, . . . Hp ^=0, 

respectively. Thus, there are altogether p(p.-l) cellB among the 
Hm H„+l 

H A complexions that contain no motes. But there is a total of p 

cells in all the possible complexions. The probability of drawing a cell 


with no motes in it is therefore 


one with one or more motes in it is 1 



and the probability of drawing 
r T 


In order to complete the solution of the problem with all kinds of 
mote present, the probability P(T,n) of drawing a oell with one or more 
motes causing freezing at a given temperature and none causing freezing 
at a higher temperature m u st be found. The probability of drawing a 
cell that would freeze at a temperature of -T° C is equal to the product 
of the probability of drawing a cell with no particles of kinds 
0,1, . . . T-l and the probability of drawing a cell with one or more 
particles of kind . T, that is, 


P(T,p) 



H 0 +Hi+H 2 + 




1 



h t 


The most probable frequency-distribution curve of spontaneous freezing 
temperature for a given value of p is' given by equation (l) . 


( 1 ) 
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Because the quantity |_i = v/AV contains both AT and V, the 
question arises as to the possible dependence of P(T,y.) on V as well 
as AV. It can be. shown that for AV/V small enough so that terms of 
the second or higher order may be neglected, the function P(T,p) is 
independent of V. Samples of volume AV ar6 drawn from two large 
volumes V and V' such that AV/v and AV/V' are very small. The 
mote densities (numbers per unit volume ) of various kinds are the same 
in both volumes, namely Pq, p^_, P 2 , . . . Pp. The corresponding expres- 
sions for P(T,p.) are 




A (Po + Pl + P2 + * * • + Pp-i)V 

1 (l «)‘* T 

7 



and 


P* = 1 


AV) 

v7 


\ (Po + Pl + P2 + * • • + Pp-l)^' 


1-1 


AV 


PrV 


With only the first two terms in the binominal expansion of 



(p 0 + p i+ p 2 + 


• + p t )v 

and 




it is found that 


P = P’ 


which proves P(T,p) to be independent of the choice of V provided 
AV/V is small enough. Even the higher-order terms are practically 
independent of V if the exponents are large compared to 1. (See the 
appendix.) 

With the expression for P(T,p) formulated in equation (l), a mote- 
density function of freezing temperature must be found such that the 
resulting probability function agrees at least in its general character- 
istics with the distribution curves presented in reference 2 . The 

+ Vi 

(T+l) tern of the geometric progression 

a, ar, ar 2 , . . . ar 1 , . . . 

where a is the number of motes in the volume V causing freezing at 
0° C and r is an arbitrary constant, was found to be a suitable jnote 
frequency-distribution function. Thus the sum of all motes causing 
freezing between 0° C and -T° C in V is 


6m 
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0 0 

The following form of equation (1) is obtained, by substitution of 
the foregoing expression: 


p( T ,„) - (ti) a 

because 


I?Q+ir^+N2+Ng+ • • 

The nature of equation (2) is evident in figure 1 where curves corres- 
ponding to particular values of p and Pq and various values of r 

are plotted. The probability curves all have a maximum and are some- 
what unsymmetrical . Furthermore, as r increases the maximum value 
of P increases and the peak of the curve becomes narrower. The curves 
of figure 1 are comparable in form to the frequency -distribution curves 
of reference 2. 

In order to fit the experimental average freezing temperature - 
droplet size curve with the theoretical curve, the mode (or 'maximum) of 
P(T,p) is found by taking the derivative dP(T,p)/dT and solving the 
equation 


% 


l-i\ i- 
r- 


1 - 


Narl 


( 2 ) 


-l' 

• +%-l = a 


o 

dT 

for T. The mode or the maximum of the probability curve, which is 
almost symmetrical, is close enough to the average so that they may be 
considered equal (table I). 

By performing the foregoing operations, the following equation for 
the mode % as a funotion of p, a, and r is obtained: 

log (log r) - log (log - log a 

T M “ log r 


( 3 ) 
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Equation (3) ia now reduced to a simple approximate form. Because V 
is much larger than AT, second- or higher-order terms in AT/T are 
neglected with the result that 


and 


JL_ 

d-1 



log Ji_ = log (1 + £V) S’ 0.4343 
fi-1 V V 


(4) 


(5) 


In the experimental curve of reference 2, Ty plotted as a function 
of the diameter D of a hemispherical drop of volume AV so that for 
comparison of theory with experiment equation (4) he comes 


JL. = _i 

“ 12V 


or 

log ^- = 0.4343 

fl _ 

Furthermore, the mote -density function is ftp = y r 1 , and the mote- 

density constant Is p 0 = . Thus, by substitution of equatlori (6) in 

equation (3) 


nD° 

12T 


( 6 ) 


T, 


M 


log (log r) - 3 log D + log 
leg r 


8.78 

P0 


(7) 


Is obtained. As a result of the assumed mote -density function, T M Is 

therefore to the first order of approximation a linear function of log D. 
The linear relation between Ty and log D is a basic one that can readily 

be compared with the experimental data of reference 2. 
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The slope of the linear equation (7) is 


T 


2 ,M~ T 1,M 


log I >2 ' - log Dj_ log r 


(a) 


where the subscripts 1 and 2 indicate two distinct values of T^ corres- 
ponding to two droplet' sizes D]_ and D 2 , respectively. The constant pq 
is found from equation ( 7 ), after setting T^ = 0, to be 


8.78 log r 

p = =-» — 

Dq 5 


( 9 ) 


where Dq is the value of D when Ty = 0. Prom equation (8), the slope 
of equation (7) is apparently a function of r alone, and Pq according 
to equation (9) determines the intercept. Equations (8) and (9) will be 
subsequently useful in evaluating the constants r and Pq from the 
slope and the intercept of the experimental average temperature -log D 
curve. Then with the constants pq and r known, the theoretical 
probability curves may be computed from equation (2) and compared with 
the experimental frequency-distribution curves of reference 2. 

In addition to equations (8) and (9) defining the mode and the inter- 
cept, an approximate transformation equation is derived in the appendix. 
With this transformation, a probability curve for one set of values of 
the parameters pq and AT can be transformed simply to a curve corres- 
ponding to another set of values. If r is kept constant and it is 
desired to transform from a distribution curve for a volume AYj_ and 


mote -density constant pq to 
given value of the probability 
Tg corresponding to AY 2 and 

to AV-l and Pq 1 "by 


another for 
P the value of the 
pq ^2 is related to 


AY2* and. Pq, 2> "then for a 
freezing temperature 
T]_ corresponding 


T 2 = Ti 


lp g (pQ,l^l/PQ,2AY 2 ) 


log r 


(10) 


This transformation is possible, because the family of curves correspond- 
ing to all values of AY and Pq are identical in size and shape but 
are displaced in the T coordinate. (See the appendix.) 
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A method of transformation from a curve for to another for 

keeping Pq and AV constant, is unavailable, but a set of curves 
for various values of r is given in figure 1 from which curves for 
other values of AY and pq may be computed using the preceding trans- 
formations. 


RESULTS AND DISCUSSION 

For comparison of the theory of this report with the experimental 
results of reference 2, the data of reference 2 for. the average freezing 
temperature -size curve are replotted on semilog paper, and an Irregular 
curve that can be-fitted reasonably well with a straight line is obtained 
(fig. 2). This result agrees with the linear relation between Tm and 
log D Indicated by equation (7). The slope of the straight line, as 
drawn in figure 2, is -7.21 and the value of Dq 1b 24 (cm). From the 
slope, the value r is found by equation (8) to be 2.60 and, from the 
intercept' Do and equation (9), Pq has the value^2 .64>&.0-4 (cm-3) . 

The linear equation that best fits the data of reference 2 is therefore 

Tjyj = -9.95-7.21 log D 

With the constants - r and P Q now known, it is possible to com- 
pute the probability curves for various droplet sizes. A comparison 
between these curves and the corresponding experimental frequency- 
distribution curves of reference 2 is given In figure 3 where histograms 
of the experimental data are reduoed to the probability scale . Fig- 
ures 3(b) to 3(d) represent the results obtained for the droplet sizes 
with the greatest amount of experimental data, and figures 3(a) and 3(e) 
represent the results for. the extremes of the size range investigated 
where the data were limited. In the computation of the probability 
curves of figure 3, a convenient value of Y = 10 3 (cm 3 ) was assumed. 

The value of Y does not Influence the curves If .it is large compared 
with AY . 


The agreement between the theoretical and experimental curves In 
figure 3 is qualitatively very good, but the following differences in 
detail exist. The modes of the theoretical curves differ from those of 
the experimental curves corresponding to the deviation of experimental 
points from the straight line drawn through them In figure 2 . The 
experimental probability curves are broader and have less pronounced 
peaks than the theoretical curves. The standard deviations (root mean 
squares of the deviations from the mean freezing temperature) of the 
theoretical curves are constant and equal to 1.38° C, but the standard 
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deviations of the experimental curves are "between 2 and 4 times as great 
and are roughly constant in the size range "below 100' microns (fig. 4). 
Above 100 microns there is a sharp inorease in the standard deviation of 

h the experimental data to a higher approximately constant value of 3.8. 

f— 1 

The foregoing difference between the experimental and theoretical 
curves cannot be ascribed to the fact that the exper ime ntal data are 
grouped into narrow droplet -size ranges to provide the freezing- 
temperature distribution curves for a given average droplet size. A 
composite theoretical curve formed by combining theoretical curves from 
the middle and the two extremes of the experimental -size group differed 
very little from the unmodified theoretical curve, as shown in fig- 
ure 3(c) . 

Better agreement between theory and experiment is obtained if 
several different mote-density functions are postulated and the corres- 
ponding theoretical curves combined. . This procedure may be justified, 
because the experimental data were taken on several days distributed 
over a period of three months. Data taken each day covered almost the 
entire range of droplet sizes investigated. If impurities were the 
active agent, the mote-density function could change from day to day 
during that time. . Thus, the combination of data taken at different 
times could easily broaden the experimental distribution curves as com- 
pared to the theoretical curves. Changes in r and pg could modify 
the slope and the Intercept, respectively, of the Ty - log D curve 
(fig. 2). 

The modes of the theoretical distribution curves may be shifted 
without changing the slope of the T^ -log D curve by varying the 
constant Pq. The weighted combination of several curves corresponding 
to values of pq between 10 "5 and 10 “3 (cm" 3 ) (see table H) results 
in theoretical curves that fit the experimental curves much better than 
the unmodified theory, as shown In figure 3. The modified theoretical 
curves, however, have a constant standard deviation equal to 2.12° C, 
which still does not agree with the standard deviations of the experi- 
mental curves, as shown In figure 4. Better agreement might have been 
obtained by varying the values of r as well as of pgj but the 
additional work involved In varying r was not considered worthwhile, 
because the variation of pg was only postulated to demonstrate the 
plausibility of the hypothesis concerning the behavior of the experi- 
mental standard deviations . The good agreement obtained between theory 
and experiment indicates that variation of the mote -density function 
among several groups of ocmbined data can account for the discrepancy 
between theory and experiment. 
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Deviations from the most probable curves may by chance occur in the 
observed distribution curves, but deviations of a given magnitude became 
less probable as the number of observations increases . This principle 
is illustrated by the greater regularity of distribution curves based on 
a large number of observed droplets as compared with those based on a small 
number of droplets. For example, the curves in figures 3(a) and 3(e) 
for -which 39 and 35 droplets, respectively, -were observed are- quite 
irregular as compared with the curve ip figure 3(c), for which 692 drop- 
lets were observed. 


CONCLUDING REMARKS 

In the derivation of the foregoing statistical’ theory, emphasis has 
been placed on the mote hypothesis of Dorsey, but the theory is also 
applicable to homogeneous nucleation as described by Volmer (reference 6) . 
According to reference 6, the random small-scale fluctuations in an 
impurity-free supercooled liquid may cause crystal nuclei to form, result- 
ing in crystal growth from the mother liquid. The presence of motes in 
the liquid would mask the homogeneous nucleation. If it were possible 
to remove all foreign nuclei, then the mote distribution function Np 
would merely be renamed the homogeneous nucleus distribution function 
and the logic of the theory in this report would apply as in the case 
of the motes except that Np would now’ be a function increasing with 
time. Thus, a container of liquid kept at a constant temperature of 
supercooling should freeze after a period of time. But- the experience 
that Dorsey (reference 1) had in maintaining three test tubes at a tem- 
perature between -8.0° and -10.3° C for 312 days might seem to invalidate 
the homogene ous -nucleation theory. The apparent inconsistency that 
exists between the results of reference 1 and those expected from the 
homogeneous -nuoleat ion theory may possibly be resolved by the following 
logic. The temperature range in which the experiment of reference 1 
was performed may be one of low nucleus -formation probability even over 
a period of 312 days. But if the probability of nucleus formation 
increases very rapidly at some temperature much lower than -10° C, the 
most probable length of time that a given water sample can remain at a 
low constant temperature of supercooling may be sufficiently reduced to 
make maintenance of the water sample in the supercooled state for a very 
long time practically impossible. 

If the spontaneous freezing temperatures of droplets are princi- 
pally determined by the presence of motes, the results of this report 
can be extended to supercooled clouds, because there should be no 
essential difference between the droplets supported on a surface and' the 
droplets suspended in the atmosphere except for possible variations In 
the mote distribution curves as a function of time and place. The 
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results of this report indicate that no well-defined temperature exists 
at which supercooled cloud droplets are transformed into ice particles. 
As the temperature is decreased in a given cloud, more and more droplets 
cm will tend to freeze until a point is reached at which very few droplets 

are left unfrozen. For example, in a cloud of uniform droplet size 
’’’’ with hemispherical diameter D of 34.5 microns or an actual spherical 

droplet diameter of 27.4 microns which is fairly common in clouds, prac- 
tically all the droplets would be frozen at a temperature of -34° C . 

On the other hand, practically all the droplets in the cloud would 
remain unfrozen down to a temperature of -18° C and only about hal f 
would be frozen at a temperature of -28° C. In clouds with a droplet 
size of 5 microns or greater, therefore, very few supercooled droplets 
exist at temperatures lower than -35° C and the existence of an icing 
cloud at such low temperatures is probably an unusual occurrence. 

The relative proportions of frozen and supercooled droplets in a 
cloud are important in determining the stability of the cloud. A cloud 
is considered to be stable if factors causing rapid dissipation of the 
cloud are absent. One of the factors that may cause a supercooled 
cloud to become unstable is the existence of ice particles in the cloud. 
Ice particles in a supercooled cloud have a lower vapor pressure than 
the surrounding supercooled drops. As a result of the vapor-pressure 
gradient existing between the liquid droplets and ice particles, water 
. tends to evaporate from the droplets and condense on the Ice particles. 
Thus, If ice particles are present in sufficient number distributed 
evenly throughout the cloud, the ice particles will grow at the expense 
of the supercooled drops until they fall out of the cloud in the form 
of precipitation (reference 7). If too few ice particles appear in a 
supercooled cloud, they will grow and fall out of the cloud leaving It 
practically intact. On the other hand, if too many ice particles appear, 
the cloud is merely transformed to an ice -particle cloud. The presence 
of ice crystals in supercooled clouds Is considered to be the principal 
factor Initiating precipitation. 

Determination of whether the statistics of freezing are affected 
by change of environment requires further Investigation. For example, 
if impurities are the cause of freezing in supercooled water, then the 
average freezing temperature of a given droplet size should be depressed 
by maintaining as .high a degree of cleanliness as possible. In addition, 
if cleanliness is found to have the anticipated influence, the effect 
of adding finely divided materials of Isjiown constitution to water can be 
investigated. In the experiments on supercooled water that have been 
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conducted up to the present time, the water has always "been in contact 
with a container or support.. In order to attain conditions that approach 
those in clouds, a method of investigating the free zing- temperature 
statistics of unsupported droplets must he found. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, May 3, 1950. 
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APPENDIX - DERIVATION OF TRANSFOIMATION EQUATION'S 

Equation (10) for transformation of a probability curve 
corresponding to one set of values of t tie parameters AY and Pq to 
another set of values is derived from equation (2) written in the form 


P(T,AV,P 0 ,r) = 


,(# 


-M 


PcpPv 


Equation (12) must first "be expanded in a series to obtain an approxi- 
mate expression, which is used to show the dependence of P on AY and 


Pq. The members of equation 


id to show the dependence of P on AY and 
(12), (l-f) P W anfi. (l - 


are written in a binomial expansion as follows! 




nT . x-nav + ^231 ■ , nr(nV-l)(nT-2 ) a£ „ /. _ 


where n represents p. 


o(^ 


or PqT t . Because' Y is very large. 


nY in general is so large that nV-1, nV-2, . . . are practically equal 
to nV so that the binomial expansion may be written as 

(l - ~ 1-nAY + ^ AY 2 - ~ AY 3 + . . . (13) 

By substitution In equation (12) of series in the form of equation (15), 

P(T,AV,p 0 ,r) ~ pqt^AY + — AY 2 - Pq 2 !^ AY 2 + 
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le obtained. When T is close to zero, r 1 is small and only the first 
term of the series in equation (14) need be used. As T increases, the 

terms of higher order must be considered and r®-l becomes more nearly 
equal to r^, because r is greater than 1. Equation (14) may there- 
fore be rewritten as 


P(T,AV,p 0 ,r) 


PqH^AY 




(15) 


The first step in the derivation of equation (10) is the proof 
that P M , the mode value of P corresponding to T M for one set 
of values of the parameters Pq and AY, is equal to % corresponding 
to another set of values. Two sets of values of p Q and AY are 
indicated by adding the subscripts 1 and 2, and the corresponding 
variables are similarly indicated. With equation (3) in the form 

log (log r) - log PrvAY - log 0.4343 

% (w. 


with values AVj and AVg of the parameter AY, with values Pq q and 
p 0,2 parameter Pq> srd- with r maintained constant, the equation 



p 0,l AV 1 


t 2,M 

r p O,2 AV 2 


(17) 


is obtained by subtraction of the equation for ^2,VL f rct!a that for 
The mode probabilities by equation (15) are 


P £ r T l,M D , 
*1,M r “ ' P 0,l 


AY X - 
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„ T 2 M (^2,11 ) 2 

’ r ’ p 0,2 4T 2 - L 


(“0,2 ^ 2 ) 2 - U^P ' (p 0 ,2 *V 2 ) 2 . . . 


Therefore, by equation (17), the separate terms in the series for ^ 
and P2,M are ©filial and 


p — p 

1,M “ 2,M 

Because for a given value of r the mode probability P^ is the 
same for all values of Pq and AV, the probability function P ranges 
from zero to P M for all values of p Q and AV. To the value of the 
probability P]_ corresponding to the values Pq AVq_, and Tj_, there 
are values Po,2j AV 2 , and T 2 such that 

P x = P 2 

From equation (15) series expansions for P-j_ and Pg are obtained anri 
substituted in this equation so that 

( ^l) 2 ( 2 

rTl p 0,l AT 1 " p 0,2 AY 2 - 2 (P0,1 A ^l) 2 " 2 (P 0 ,2 AT 2) 2 + 


P 0 ,l AY l ) 2 - 


(P 0,2 AY 2) 2 


For values of T close to zero, the probabilities and P 2 are given 

to a good approximation by the first term of the series in equation (15), 
and the right-hand member of equation (18) is therefore negligible. As T 
increases, the terms of second and higher older in equation (15) become 
appreciable, and the individual terms on the right-hand side of equa- 
tion (18) are no longer negligible. But as T approaches T M , the 

bracketed differences on the right-hand side of equation (18) approach 
zero according to equation (17). The differences on the right-hand side 
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of equation (18) are therefore negligible in the range of T from. 0 
to a value somewhat greater than T^. Therefore, because P becomes 

negligible for T Jus t“a- little greater than Tj^ (see fig. 3), the 
equation 


T-i ^ Tp 

r 1 p AV = r 2 

M 0,l 1 


p AY 

0,2 2 


(19) 


is a good approximation in the range of T where P is of-appreciable 
magnitude. If the log of equation. (19) is taken, this equation may be 
rewritten as a transformation equation, which yields the displacement 
of a probability curve in the T coordinate as 


T 


2 


log (p Q>1 Av l/ p Q ,2 AY 2 ^ 

T 1 + - log r 


( 10 ) 


Calculations of P for various values of p Q and AY from equa- 
tion (2) confirm the approximate validity of equation (10). 
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TABLE I - COMPARISON OF MODES AND AVERAGES OF 
FREQUENCY -DISTRIBUTION CURVES 


Diameter 
of hemi x 
spherical 
drops 
(microns ) 

Mode 

(°c) 

Average 

(°c) 

(Mode )- (average ) 

(°c) 

llxperimental' values 

805 

-16.6 

-17.4 

0.8 

161 

-23.3 

-24.2 

.9 

69 

-25.5 

-25.2 

-.3 

34.5 ‘ 

-24.4 

-26.2 

1.8 

8.75 

-32.2 

-32.1 

-.1 

Theoretical values 

1564 ■ 

-15.8 

-15.2 

-0.6 

156.4 

-23.0 

-22.6 

-.4 

72.6 

-25.5 

-24.9 

-.6 

33.70 

-27.7 

-27.3 

-.4 

7.26 

-32.7 

-32.1 

-.6 


lata from, reference 2. 


TABLE U - VALUES OF p Q AND CORRESPONDING WEIGHTS 
USED IN MAKING COMPOSITE CURVE 


P 0 

Weight 

(cm - ^) 


1.0 x io" 3 

2 

2.64 X 10~ 4 

4 

1.05 X 10“ 4 

2 

5.7 X 10 -5 

2 

3.7 x 10" 5 

2 

2.6 X 10 -5 

1 
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Figure 1. - Probability curves for various values of r < 
1-AV/V = (n-l)/n * 1 - 10" 6 , AV - 10 -3 , V * 10 J 
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-10 ^14 -18 - -22 -26 

Freezing temperature, °C 


I I I 1 

10 0 . . -1C . . . .-20 

Freezing temperature, °F 

(a) Diameter, D, 805 microns. 

Figure 3. - Comparison of theoretical and experimental probability ourvea. 
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10 0 -10 ■ ■ -20 -50 

Freezing temperature, °F 

(o) Diameter, D, 69 microns. 

Figure 3. - Continued. Comparison of theoretical and experimental 

probability curves. 
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Freezing temperature, °F 
(d) Diameter, D, 34.5 microns 

Figure 3. - Continued^ Comparison of theoretical and experimental 

.probability curves. 
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I I I I I I 

10 .0 ' -10 ... ... -80 _. -30 -40 

■ Freezing temperature ; °F 

(e) Diameter, D, 8.75 microns. 


Figure 3. - Concluded. Comparison of theoretical and experimental 

probability curves. 
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Figure 4, -Comparison of standard devlationa of theoretical distribution curves with thoae of 

experimental curves. 
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